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We report polarized- and unpolarized- neutron inelastic scattering measurements of the magnetic 
excitation spectrum in the spin-charge ordered phase of La 3/ / 2 Sr 1 / 2 Ni04. Up to energies of ~30 meV 
we observe broad magnetic modes characteristic of a near checkerboard ordering. A linear spin- 
wave model for an ideal checkerboard ordering with a single antiferromagnetic exchange interaction 
J' = 5.8 ± 0.5 meV between next-nearest-neighbour spins on Ni 2+ sites, together with a small 
J\ F-like single-ion anisotropy, provides a reasonable description of the measured dispersion. Above 
30meV the excitations are not fully consistent with the linear spin- wave model, with modes near 
the two-dimensional reciprocal space wavevector (0.5, 0.5) having an anomalously large intensity. 
Furthermore, two additional dispersive modes not predicted by spin wave theory were observed, both 
of which are probably magnetic. One disperses away from (0.5, 0.5) in the energy range between 
50-56 meV, and the other appears around (h, k) type positions (h, k = integer) in the energy range 
31-39 meV. We propose a model in which these anomalous features are explained by the existence 
of discommensurations in the checkerboard ordering. At low energies there is additional diffuse 
scattering centred on the magnetic ordering wavevector. We associate this diffuse scattering with 
dynamic antiferromagnetic correlations between spins attached to the doped holes. 



I. INTRODUCTION 

The existence of charge ordering phenomena continues 
to attract considerable interest, especially concerning the 
properties of layered transition- metal oxides. The na- 
ture of the charge ordering in these materials depends on 
the degree of doping relative to a Mott insulating phase. 
At low doping the charges tend to segregate into paral- 
lel lines. The resulting stripe superstructures have been 
widely discussed in connection with high-temperature su- 
perconductivity in cupratesi At higher doping there is a 
competition between charge-ordered and metallic states. 
The existence of magnetic correlations is also an impor- 
tant factor, and the coupling between charge, orbital, 
spin and lattice degrees of freedom can often stabilize 
novel ordered states. 

At half doping, systems with dominant Coulomb re- 
pulsions are expected to exhibit a stable Wigner crystal 
state, which in the case of a two-dimensional square lat- 
tice takes the particularly simple form of a checkerboard 
pattern. This type of charge order has been observed 
in the isostructural layered perovskites La 3 / 2 Sr 1 / 2 Ni04, 2 



La 3 / 2 Sr 1 / 2 MnO 



4^ 



and La 3 / 2 Sr;iy 2 CoO 



Interest- 



ingly, checkerboard charge order has also recently 
been observed at a doping levels well below 0.5 in 
Bi 2 Sr 2 CaCu20 8+ 5^, and Ca 2 _ a: Na a ;Cu0 2 Cl 2f & with a 
periodicity of 4 Cu atoms. Further experimental stud- 
ies of simple systems with checkerboard charge order 
could help explain the formation of electronically-ordered 
states in more complex systems. 

In this paper we describe a neutron inelastic scat- 



tering study of the spin excitations in the half-doped 
layered nickelate La 3 / 2 Sr 1 / 2 Ni04. This compound ex- 
hibits checkerboard charge order below Tco ~ 480 K. 2,7 
Below Tic ~ 180 K the checkerboard pattern becomes 
slightly incommensurate, and below T$o ~ 80 K it is ac- 
companied by incommensurate magnetic orderi A spin- 
reorientation transition has been observed at Tsr, = 
57 K. 8 By mapping out the spin excitation spectrum of 
La 3 / 2 Sr 1 / 2 Ni04 as a function of energy and wavevector 
and comparing the dynamic response with that of other 
charge-ordered compounds we hope to gain a better un- 
derstanding of the interactions that stabilize the ground 
state in these systems. 

Previously, the spin excitation spectrum of 
La 2 _ a; Sr a ;Ni04 has been investigated over a wide 
range of wavevector and energy in stripe-ordered com- 
pounds with x = 0.275,2 x = 0.31^ and x = l/3iii±2 
The spin excitations are highly two-dimensional (2D), 
and fairly independent of x over this small range of 
compositions studied. Over most of the energy range 
the magnetic dynamics can be described in terms of spin 
wave excitations of the antiferromagnetically ordered 
Ni spins in the region between the charge stripes^i At 
low energies, however, a second magnetic component is 
observed that has been shown to come from quasi- ID 
antiferromagnetic spin fluctuations along the stripes 
themselves 



Figure ^ illustrates how the different ordered phases 
in La 3 / 2 Sr!/ 2 Ni04 can be identified in diffraction mea- 
surements. Figure da) is a simplified model of the 
ground state spin-charge order within the Ni0 2 layers 



of La3/2Sr!/ 2 Ni04 neglecting the small incommensurate 
modulation. For present purposes we assume the doped 
holes are site-centred, so that the checkerboard pattern 
is derived from alternating Ni 2+ and Ni 3+ ions carrying 
spins S = 1 and S = 5, respectively. The S — 1 spins 
are assumed to be ordered antiferromagnetically, but we 
make no assumption at this stage about magnetic order 
of the S — h spins (this issue will be addressed later in 
the paper). The positions in reciprocal space of the cor- 
responding Bragg peaks are shown in Fig. ^b). Peaks 
from the charge order have two-dimensional wave vec- 
tors (in units of 2k /a) (h + \,k + i), where h and k 
are integers^ The magnetic order has double the pe- 
riodicity of the charge order, so peaks from the mag- 
netic order appear at (h + 5, k + 5) ± (|, |). Rotation 
of the ordering pattern by 90deg generates an equiva- 
lent magnetic structure, this time with magnetic peaks 
at fc+|)±(i — j). In the absence of a symmetry- 

breaking interaction we expect an equal population of 
these two domains, so the pattern of Bragg peaks will be 
a superposition, as shown in Fig. ^b). 

The actual spin-charge ordered phase of 
La 3 / 2 Sr 1 / 2 Ni04 observed below Tic does not con- 
form to the ideal structure shown in Fig.^a). Instead, 
the magnetic Bragg peaks are found at the incommen- 
surate positions (h + h, k + |, I) ± (e/2, e/2, 0) with I an 
odd integer, and (h + |, k + |, I) ± (e/2, -e/2, 0) with I 
an even integer, where e « 0.44. 7-8 -i£ New charge-order 
satellite peaks appear at (h ± e, k ± e) in addition to the 
checkerboard charge order peak at (h + 4, k + h), with 
little or no I dependence*^ The full set of 2D magnetic 
and charge order wavevectors for the incommensurate 
phase of Las^Sr^NiC^, including those for the 90 deg 
domain, is shown in Fig. ^c). 

Kajimoto et al. 7 suggested that two types of charge 
order coexist in the incommensurate phase, one part of 
the system being charge-ordered in a checkerboard pat- 
tern and the other part adopting an incommensurate, 
stripe-like, spin-charge order. These authors developed 
models for the latter component based on the introduc- 
tion of diagonal discommensurations in the ideal checker- 
board spin-charge structure of Fig. ^a). As pointed out 
by Kajimoto et al., the stability of the incommensurate 
structure is probably the result of a competition between 
magnetic and electrostatic energy. The strong superex- 
change interaction favours having antiparallel spins on 
nearest neighbour Ni sites, whereas the Coulomb inter- 
action tries to have a uniform charge density. 



II. EXPERIMENTAL DETAILS 

The neutron scattering measurements were performed 
on a single crystal of La3/2Sr 1/ / 2 Ni04 grown by the 
floating-zone method^ The crystal was a cylinder 35 mm 
in length and 7 mm in diameter. The oxygen content of 
the crystal was determined by thermogravimetric analy- 
sis (TGA) to be 4.02±0.01. This is the same crystal that 
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FIG. 1: (Color online) (a) Ideal checkerboard spin-charge or- 
dering in a Ni02 square lattice. Circles represent Ni 3+ holes, 
and solid arrows represent spins on Ni 2+ sites. The broken 
lines are included to highlight that the spin pattern breaks 
the 2D symmetry of the checkerboard charge ordering. J' 
is the exchange coupling parameter of the Ni 2+ spins across 
the Ni 3+ site. This commensurate ordering is not realized 
in practise in La 3/ / 2 Sr 1/ / 2 Ni04. (b) Diagram of part of the 
(h, k) plane in 2D reciprocal space showing the positions of 
the first-order magnetic and charge order Bragg peaks for the 
ideal checkerboard ordering represented in (a). The peaks 
from the equivalent domain in which the magnetic ordering 
is rotated by 90 deg relative to that in (a) are superimposed, 
(c) The same diagram as (b) except with the magnetic and 
charge order Bragg peaks observed in the incommensurate or- 
dered phase of Laa^Sr^NiC^. For simplicity we neglect the 
variation in the peak positions in the direction perpendicular 
to the Ni0 2 plane. 



was used in our earlier neutron diffraction and magneti- 
zation study of the magnetic order£ 

The majority of the unpolarized-neutron scattering 
measurements were performed on the time-of-fiight chop- 
per spectrometer MAPS at the ISIS Facility. The crystal 
was mounted on MAPS in a closed-cycle refrigerator and 
aligned with the c axis parallel to the incident beam di- 
rection. A Fermi chopper was used to select the incident 
neutron energy. Incident energies of 60 and 100 meV were 
used. The intensity was normalized and converted to 
units of scattering cross-section (mbsr -1 meV -1 [f.u.] -1 ) 
by comparison with measurements from a standard vana- 
dium sample. Scattered neutrons were recorded in large 
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banks of position-sensitive detectors. The spin disper- 
sion was found to be highly two dimensional. Hence, we 
analyzed the data by making a series of constant energy 
slices and projecting the intensities onto the (h, k) two- 
dimensional reciprocal lattice plane. 

Further unpolarized- and polarized-neutron measure- 
ments were performed on the triple-axis spectrometers 
IN8 and IN20 at the Institut Laue-Langevin. The ener- 
gies of the incident and scattered neutrons were selected 
by Bragg reflection from a double-focusing bent Si crys- 
tal monochromator on IN8, and an array of Heusler alloy 
crystals on IN20. On both instruments the data were 
obtained with a final neutron wavevector of 2.66 A -1 , 
and a pyrolytic graphite (PG) filter was placed between 
the sample and analyzer to suppress higher-order har- 
monic scattering. For polarized-neutron scattering on 
IN20 the spin polarization, P, was maintained in a spec- 
ified orientation with respect to the neutron wavevector, 
Q, by an adjustable guide field of a few mT at the sam- 
ple position. For the experiment on 1N8 the crystal was 
orientated with the [100] and [010] crystal directions in 
the horizontal scattering plane, so that (h, k, 0) posi- 
tions in reciprocal space could be accessed. On IN20, we 
mounted the crystal with the [001] and [110] directions in 
the horizontal scattering plane, so that (h, h, I) positions 
in reciprocal space could be accessed. 



III. RESULTS 
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Figures[2Ia)— (c) show constant-energy slices from runs 
performed at T = 10 K on MAPS. 1 - The slices have 
been averaged over a range of energies, as indicated in 
the figures. In (a), the range is 4-6 meV, and the in- 
tensity is seen to be enhanced at the four equivalent 
magnetic order wavevectors (0.5,0.5) ± (e/2, e/2) and 
(0.5,0.5) ± (e/2, -e/2), where e ps 0.44. This signal, 
therefore, corresponds to low energy spin excitations 
from the magnetically ordered ground state. As the en- 
ergy increases, the intensity spreads out but initially re- 
mains centred on the magnetic wavevectors. This can 
be seen in Fig.[3[b), which shows data averaged over 25- 
30 meV. At still higher energies the four blobs of intensity 
tend to merge towards (0.5, 0.5). This is illustrated in (c), 
which corresponds to an energy range 40-45 meV. 

In Figs. Efd)— (f) we plot cuts through the data in the 
(£,£) direction for the same energy ranges as used in 
Figs. |5Ja)-(c). At low energies, E = 4 — 6 meV, the 
scattering takes the form of a single peak centred on the 
magnetic wavevector. The fitted Gaussian half- width of 
this peak converts to a correlation length of 35 ± 1 A(the 
inverse of the half width at half maximum). At energies 
of 25-30 meV [Fig.|2Je)] the lincshape is also centred on 
the magnetic wavevector, but it now contains two peaks 
which can just be resolved [another example can be seen 
in Fig.[5Jg)]. The fitted width of these peaks corresponds 
to a correlation length of 24±2 A. As the energy increases, 
the right-hand of the two peaks grows while the left-hand 



FIG. 2: (Color online) Neutron scattering measurements from 
La 3/2 Sr 1/2 Ni04 measured on MAPS at T = 10 K. (a) - (c) 
are constant-energy slices showing the variation of the inten- 
sity in the (h, k) plane at different energies. The data are 
averaged over the range of energies indicated above the fig- 
ures. The incident neutron energy was 60meV for (a), and 
100 meV for (b) and (c). Data from four equivalent Brillouin 
zones have been averaged, (d) — (/) show cuts along the (£, £) 
direction for the same energies as shown in (a) — (c). The solid 
lines are the results of fits with one Gaussian [(d) and (f)] or 
two Gaussians [(e)] on a linear background. In (d) we also 
display the (0.275, 0.275, 5) magnetic order Bragg peak mea- 
sured on IN20. The out of plane wavevector for the position 
(0.275, 0.275) are as follows (o)&(d) is I = 0.55, (6) I = 2.3, 
(e) I = 3.0, (c) I = 3.45 and (/) I = 5.2. 



peaks diminishes. Above E w 40 meV the left-hand peak 
has virtually no intensity, as can be seen in Fig. Elf). 

Figure |21 shows the dispersion of the spin excitations 
obtained from Gaussian fits to the peaks observed in cuts 
such as those displayed in Fig. Etd)-(f). The points in 
Fig. |21 are the fitted peak centres corrected for the effect 
of the non-zero width of the cut perpendicular to the cut 
direction. For this purpose the scattering was taken to 
be a circle centred on the magnetic wavevector. Within 
the experimental limits of our data the dispersion was 
found to be the same in orthogonal cut directions. 
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FIG. 3: (Color online) Dispersion of the magnetic excitations 
in La3/2Sr!/ 2 Ni04 parallel to the (£, £) direction. The points 
are the results of fits to cuts such as those shown in Fig.|3 The 
broken line is calculated from Eq. JSJ, the spin wave dispersion 
for a checkerboard ordered system with exchange parameter 
J' = 5.8 me V and out-of-plane anisotropy K c — 0.05 me V. 
Triangles above the calculated dispersion curve indicate the 
positions of an additional observed scattering mode. The 
shaded area represents the region in which additional scat- 
tering hampers the study of the spin-wave excitations, see 
Fig.il 
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(Color online) Neutron scattering measurements 
. 3/2 Sr 1/2 Ni04 measured on MAPS at T = 10 K. (a) 



FIG. 4 

from La ; 

Constant-energy slice through the data showing the intensity 
distribution in the (h, k) plane in the energy range 51-54 meV. 
The slice reveals weak diffuse scattering symmetrically dis- 
tributed around (0.5,0.5). Data from four equivalent Bril- 
louin zones have been averaged to improve the statistics, (b) 
Cut through the 51-54 meV slice along the path indicated by 
the arrow in (a). 



In constructing the dispersion curve we had to take 
care to avoid confusing the scattering from magnetic ex- 



citations and phonons, especially in the energy range 
10-40 meV where the phonon scattering is particularly 
strong. One check we made was to compare data ob- 
tained at 300 K with that obtained at 10 K. Phonon scat- 
tering increases in strength with temperature, whereas 
magnetic scattering decreases in strength above the mag- 
netic ordering temperature. 

At energies above 50 meV we observed a broad ring of 
scattering centred on (0.5, 0.5). This is illustrated in Fig. 

which shows (a) a constant-energy slice averaged over 
the energy range 51-54 meV, and (b) a cut through this 
data in the (£, 1 — direction. This mode disperses away 
from (0.5, 0.5) with increasing energy, but was too weak 
to measure above 56meV. This scattering was found to 
be slightly weaker at 300 K than at 10 K, 34± 19 origin for 
it. We have added the peak positions of this scattering 
to the dispersion curve in Fig. |3] 

In the energy range 31 — 39meV we observed an in- 
teresting scattering feature separate from the main spin 
wave scattering. Figure |5f a) shows the distribution of 
intensity measured on MAPS for energies averaged over 
the range 33-37 meV covering a large area of reciprocal 
space. This intensity map displays strong diffuse scat- 
tering rings centred on the reciprocal lattice zone cen- 
tres (h,k), where h and k are integers. Figures EJb)- 
(d) present constant-energy slices through this scatter- 
ing centred on energies of 32meV, 34meV and 36meV. 
In these slices the symmetry-equivalent data have been 
folded into one quadrant to improve the statistics. The 
intensity is seen to disperse away from the (1,0) point. 
In these measurements the out-of-plane component of 
the scattering vector varies from I = 2.5 to 4. A search 
was made with the IN8 spectrometer of other Brillouin 
zones in the (h, k, 0) reciprocal lattice plane. This survey 
confirmed the results found on MAPS, showing that the 
scattering was not restricted to a particular out-of-plane 
component of the scattering vector. 

To assess whether or not this scattering is magnetic 
we performed measurements as a function of tempera- 
ture. Cuts through the MAPS data at 10 K and 300 K 
are shown in Fig. |5{e) . These data show the scattering to 
decrease slightly with temperature, being approximately 
45 ± 12% stronger at 10 K than at 300 K. We also ob- 
served similar scattering in the range 25-30 meV. Fig|Sff) 
shows Q scans through the point (2, 0, 0) at an energy of 
25meV for temperatures of 1.5 K and 300 K, measured 
on IN8. The scattering here appears to be temperature 
independent. To understand these temperature effects 
we need to take into account that this diffuse scattering 
ring overlaps the magnetic ordering wavevectors, so will 
contain some spin-wave scattering. The spin wave scat- 
tering itself will, of course, be temperature dependent. 
To illustrate this we show in Fig. Elg) scans through 
the spin-wave scattering associated with the magnetic 
wavevector (1.5,1.5) — (e/2,e/2) measured at a slightly 
higher energy (E = 30meV) on IN8 at T = 10 K and 
300 K. The amplitude of the spin wave scattering clearly 
decreases with temperature. From this loss of intensity 
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FIG. 5: (Color online) Excitations from La3/2Sri/2Ni04. (a) 
Distribution of scattering intensity in the (h, k) plane mea- 
sured on MAPS at T = 10 K. The scattering has been aver- 
aged over the energy range E — 33 — 37meV. (b)-(d) Slices 
through the scattering intensity measured at T = 10 K in the 
(h , k) plane averaged over the energy ranges (b) 31 — 33 meV, 
(c) 33 - 35 meV, (d) 35 - 37 meV. The data in the MAPS de- 
tector have been folded to take advantage of the symmetry of 
the scattering, (e) Cuts through the data shown in (c) along 
a line parallel to the (h, h) direction passing through (1, 0). 
The data at T = 10 K and 300 K, have been offset by 5 units. 
The lines in (e) are the result of fitting the points with two 
gaussian functions on a sloping background, (f) Scans of the 
scattering observed at E = 35 meV for T = 10 K and 300 K 
taken on IN8. (g) A scan through the spin wave dispersion 
at E = 30 meV for T = 10 K and 300 K taken on IN8. The 
300 K data has been offset by the addition of 1000 counts. 



we can infer the temperature dependence of the remain- 



ing component of the scattering in Figs. E{a) -(f). Our 
analysis indicates that at E ~ 35meV the non-spin- 
wave component of the diffuse scattering ring presented 
in Figs.|3Ja)-(d) decreases in intensity with temperature, 
whereas at E = 25meV [Fig.[5[f)] it increases with tem- 
perature. This leads us to conclude that the diffuse scat- 
tering around the Brillouin zone centres shown in Figs. 
EJa)— (d) is most likely magnetic in origin. 

Further support for this conclusion was provided 
by examining data we have collected on crystals of 
La2- 2; Sr 2 ;Ni04 with x — 1/3 and 0.275 under simi- 
lar condition during a separate experiment on MAPS. 18 
Constant-energy slices for x — 1/3 and 0.275 in the en- 
ergy range 31 — 39 meV do not show any diffuse scattering 
rings like those shown in Figs. OJa)-(d), so this feature 
seems to be specific to x = 1/2. This makes it highly un- 
likely that this scattering comes from phonons intrinsic 
to the host lattice of I^-aSr-cNiO^. 

For energies below 30 meV we performed additional 
measurements with polarized neutrons. In this energy 
range there is strong scattering from phonons, and po- 
larized neutrons were necessary to provide an unambigu- 
ous separation of the magnetic and non-magnetic scat- 
tering. We were particularly interested in studying the 
magnetic scattering as a function of energy because ear- 
lier measurements on I^-zSr^NiC^ compositions with 
x = 0.275, 0.33 and 0.37 had revealed unexpected struc- 
ture in the energy spectrum^iiS For these measurements 
the neutron polarization P was aligned parallel to the 
scattering vector Q so that the spin-flip (SF) channel 
would contain purely magnetic scattering. 

The most interesting finding is reproduced in Fig. 
which shows energy scans measured in the SF channel on 
samples with compositions x — 1/3 and x = 1/2 at their 
respective magnetic ordering wavevectors. The x = 1/3 
data contains two peaks, one centred on 7meV and the 
other on 26 meV. The x = 1/2 data contains a single peak 
at 5 meV. The lower energy peaks were shown to be gaps 
due to spin anisotropy i 11 ' 19 Below these gaps the out-of- 
plane spin fluctuations are quenched. The origin of the 
higher energy peak in the x = 1/3 data, which was first 
discussed in Ref. ^] and which was also found in sim- 
ilar data from crystals with x = 0.275 and x = 0.37, 19 
remains a mystery. Of all the compositions studied, the 
one with x — 1/3 exhibits this feature most strongly. 
The absence of a corresponding peak in the x = 1/2 
data seems to suggest that this peak is a property of 
compounds with static stripe order with a periodicity of 
^3 lattice spacings. The monotonic decrease in intensity 
above 5meV in the x = 1/2 data is qualitatively con- 
sistent with the expected 1/E dependence of the cross 
section for scattering from antifcrromagnetic spin waves. 

Finally, we describe a feature observed in the scatter- 
ing from La3/ 2 Sr 1 / 2 Ni04 at low energies. Figure |JJ a) is 
a map of the intensity measured on IN8 covering part 
of the (h, k, 0) plane in reciprocal space. The map was 
constructed from a series of scans performed parallel to 
(h, h, 0) at an energy of 3 meV and a temperature of 2 K. 
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FIG. 6: (Color online) Polarized-neutron scattering from 
charge-ordered L^-^Sr^NiCU with x = 1/3 and 1/2. The 
plots show energy scans performed at the magnetic ordering 
wavevector appropriate to each composition. The spin-flip 
(SF) scattering was measured with P || Q, and the SF back- 
ground was estimated from scans centred away from the mag- 
netic ordering wavevector. The measurements were made at 
10 K for x = 1/2 and 13.5 K for x = 1/3. 



Strong scattering can be seen centred on the magnetic 
ordering wavevectors, but additional weak diffuse scat- 
tering can also be seen centred on these same positions. 
This diffuse scattering is slightly elongated in the diag- 
onal directions, parallel to the discommensuration lines 
in the distorted checkerboard structure. There was no 
observable elastic diffuse scattering, so this feature rep- 
resents a short-range dynamic magnetic correlation. We 
were able to follow the diffuse inelastic scattering up in 
energy to ~10meV. 

In Fig. E^b) and (c) are shown scans through the spin 
wave scattering and the diffuse scattering made along the 



lines marked A and B, respectively, in Fig. Ufa). Figure 
E{d) displays the temperature dependence of the inte- 
grated intensity of the peaks in these two scans. The 
intensity for scan A is seen to increase with increasing 
temperature. This is due partly to thermal population 
of the spin waves, and partly to the reorientation of the 
ordered moment of the Ni 2+ spins that takes place at 
57 K£ Scan B, on the other hand, shows an initial in- 
crease in intensity on warming which peaks at around 
20 K before decreasing at higher temperatures. This de- 
crease suggests that the diffuse scattering is magnetic in 
origin, and the striking difference in temperature depen- 
dence between the spin wave scattering and the diffuse 
scattering strongly indicates that they arise from two dif- 
ferent magnetic components in the system. 



IV. DISCUSSION 

Let us first summarize the key observations, and then 
provide some quantitative analysis. From the various 
neutron scattering measurements we have performed on 
spin- and charge-ordered La3/ 2 Sr 1 / 2 Ni04 we have been 
able to identify four distinct features: (1) dispersive 
spin excitation associated with the magnetic ordering 
wavevector; (2) low energy diffuse spin dynamics also 
associated with the magnetic ordering wavevector but 
with a distinct temperature dependence; (3) a proba- 
ble magnetic mode dispersing from (0.5, 0.5) at ener- 
gies in the range 50-56 meV; (4) a probable magnetic 
mode dispersing from (h, k) type positions in the energy 
range 31-39 meV. All these features are relatively broad 
in wavevector, and therefore arise from dynamic correla- 
tions that are short-range in nature. 

An obvious starting point for any analysis of the mag- 
netic excitations in I^^Sr^NiC^ is the ideal checker- 
board spin-charge ordering pattern shown in Fig. ^a) . 
As far as the spins are concerned, this is a simple two- 
sublattice antiferromagnet built from 5 = 1 spins on the 
Ni 2+ sites. Spins attached to the Ni 3+ sites are ignored 
for the time being. Following our previous work on stripe- 



ordered La 5 / 3 Sr 1 / 3 NiO 
of the form 



li 



4. 



we adopt a spin Hamiltonian 



H = J' 

(id) 



z\2 



(1) 



where the first summation describes the exchange in- 
teractions between pairs of Ni 2+ spins in linear Ni 2+ - 
0-Ni 3+ -0-Ni 2+ bonds, and the second summation de- 
scribes the small, AF-like, single-ion anisotropy. Here, 
as in Ref. [ill J' is defined as the exchange energy per 
spin (multiply by 2 to obtain the exchange energy per 
bond). We neglect the diagonal exchange couplings be- 
tween Ni 2+ sites which are needed to stabilize the spin 
arrangement but are assumed to be small relative to J'. 
In effect, therefore, the system is treated as 2 uncou- 
pled square-lattice antiferromagnets with lattice param- 
eter 2a. 
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FIG. 7: (Color online) Low energy scattering from 
La3/ 2 Sr 1 / 2 Ni04 measured on IN8 at an energy of 3 meV and a 
temperature of 2 K. (a) Map of the scattering intensity in the 
(h, k, 0) reciprocal lattice plane, (b) and (c) Cuts along the 
lines marked A and B in (a), (d) Temperature dependence 
of the integrated intensity of the peaks shown in (b) and (c) . 
The lines are guides to the eye. 



The magnon dispersion derived from Eq. JJ| in the 
linear approximation is given by 



7(Q) = -j^osi^Qxa) +cos(2Qj,a)]. 



(3) 



E(Q) = 8J'S{(l+K c /8J / ) 2 ~h(Q)±K c /8J'} 2 } 



7/12T.1/2 



(2) 



The splitting of the two branches of the dispersion curve 
is such that at the magnetic zone centre one mode is 
gapped and the other isn't. The size of the gap is 
4S{2J'K C ) 1 ' 2 . When K c < J' the maximum energy of 
the dispersion curve is approximately 8 J'S. 

In Fig.[3]we have plotted the spin wave dispersion along 
the direction calculated from Eq. with S = 1, 

J' = 5.8 meV and K c — 0.05 meV. These parameters 
were chosen to match the observed maximum spin wave 
energy (~ 45meV) and anisotropy gap (~ 3meV). The 
spin- wave dispersion curve is seen to provide a reasonable 
description of the experimental data, apart from the ob- 
vious shift from the observed incommensurate wavevec- 
tor (0.275, 0.275) to the ideal checkerboard wavevector 
of (0.25,0.25). There is no detectable scattering from 
spin wave modes on the dispersion curve near £ = and 
£ = 0.5 because the antiferromagnetic structure factor is 
small in the magnetic Brillouin zones centred on (0, 0) 
and (0.5,0.5). 

From this analysis we can give a rough estimate of the 
exchange and anisotropy parameters for La3/ 2 Sr 1 / 2 Ni04. 
After consideration of the experimental errors, these are 
J' = 5.8 ± 0.5 meV and K c = 0.05 ± 0.02 meV. It 
is interesting to compare these values with those de- 
rived from similar spin- wave analyses performed on other 
La2_ a; Sr a ;Ni04 compositions. For x = 1/3 the exchange 
parameters were found to be J = 15 ± 1.5 meV, J' = 
7.5 ± 1.5 meV and K c = 0.07 ± O.OlmeVfii where J is 
the exchange interaction between Ni 2+ spins on nearest- 
neighbour lattice sites. For undoped La2Ni04 the results 
were J = 15.5 meV and K c = 0.52meV, 2Q This com- 
parison shows that J' and K c are similar for x = 1/3 
and x = 1/2, but that K c is very much larger in un- 
doped La2Ni04. An explanation for why the single-ion 
anisotropy reduces so dramatically with doping is so far 
lacking. 

We now discuss some of the obvious shortcomings of 
the model. We have already mentioned that the mag- 
netic ordering wavevector in La3/2Sri/2Ni04 is q m = 
(0.275,0.275) not (0.25,0.25). Other problems are: (1) 
the spin-wave scattering intensity above ~ 25 meV be- 
comes progressively more asymmetric about q m with 
increasing energy, i.e. stronger on the side nearest to 
(0.5, 0.5) — see Figs.EIe) and (f), and Fig.^g). This dis- 
agrees with the spin wave theory for the model described 
above, which predicts a symmetric scattering intensity 
about the magnetic zone centre. (2) The extra scatter- 
ing intensity observed around (0.5,0.5) above 50meV, 
and around (1,0) and equivalent positions in the energy 
range 31 — 39meV suggest extra magnetic modes not 
present for a checkerboard ordering. (3) The spin wave 
scattering is very broad, implying that the spin waves 
propagate only a few lattice spacings before being scat- 
tered or decaying into another excitation channel. (4) 



The source of the low energy diffuse scattering needs to 
be identified. 

To explain these features unequivocally we need a com- 
plete description of the static ordering, which we don't 
have. However, the idea that the static incommensurate 
order in L^^Sr^NiC^ can be understood in terms of 
an ideal checkerboard pattern broken up periodically by 
discommensurations is physically appealing^ and might 
provide some clues as to the origin of the various features 
in the spin excitation spectrum. 

The two simplest types of discommensuration in the 
checkerboard pattern are illustrated in Fig. |S] The first 
[Fig. IBfa)] has a line of nearest-neighbour antiparallel 
spin pairs coupled by a superexchange interaction J ex- 
pected to be close to 15meV as found in the x = 
and x = 1/3 compounds. Around the discommensura- 
tion the local hole density is below average. The second 
[Fig. OJb)] has a line of holes on oxygen sites, which in- 
creases the local hole density and forces parallel align- 
ment on neighbouring Ni spins through the double ex- 
change mecahanism. We will refer to these as antiferro- 
magnetic (AFM) and ferromagnetic (FM) discommensu- 
rations. More complex discommensurations involving a 
greater degree of perturbation of the ideal checkerboard 
pattern are also possible. 

Discommensurations provide a mechanism for produc- 
ing an incommensurate spin density wave modulation of 
the checkerboard pattern, as discussed earlier. At low 
energies, the magnon dispersion is expected to be sim- 
ilar to that of the ideal checkerboard antiferromagnetic 
ordering, except that the magnon dispersion is shifted 
away from (0.25,0.25) to the incommensurate wavevec- 
tor, as observed. 

Let us focus, however, on the specific structure of the 
discommensurations themselves, which is expected to in- 
fluence the spin excitation spectrum at higher energies. 
The AFM discommensuration can be regarded as a zig- 
zag chain with AFM intra-chain exchange J coupled to 
the checkerboard AFM background. The latter has ex- 
change J', and J w 2J'. Above the maximum energy of 
the spin wave dispersion (~ 45 meV) we might expect to 
observe spin excitations characteristic of an AFM zig-zag 
chain. Because there will be discommensurations running 
along both diagonals the scattering will take the form of 
a square of scattering dispersing away from (0.5, 0.5) type 
positions, where the structure factor of the AFM zig-zag 
chain is largest. This could explain the observed 'ring' of 
scattering apparently dispersing from (0.5,0.5) at ener- 
gies above 50 meV. At lower energies the spin excitations 
will have mixed checkerboard and zig-zag chain charac- 
ter, and this could explain why the spin wave scattering 
becomes stronger on the side nearest to (0.5, 0.5). 

Similarly, FM discommensurations resemble FM zig- 
zag chains. Spin excitations with FM zig-zag chain char- 
acter are expected to be observable near to the structural 
zone centres. The strength of the intra-chain double ex- 
change is not known, but this effect could account for 
the mode dispersing from [h, k) type positions observed 
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(b) FM Discommensuration 



Oxygen Centred 
Stripe 



b 



FIG. 8: (Color online) Two types of discommensuration that 
could exist in Las^Sr^NiO^i- Circles and arrows denote 
holes and 5* = 1 spins, respectively, on the Ni sites, in (a) the 
couplings between next nieghbour Ni 2+ spins, J, and Ni 2+ 
spins across the Ni 3+ sites, J', are indicated. The circled x 
symbols show the positions of an oxygen centred charge stripe. 
The spins on nearest neighbour Ni sites have (a) antiparallel 
and (b) parallel alignment, respectively, in the discommen- 
suration (indicated by a dashed line). In (b) holes occupy 
oxygen sites in the discommensuration. 



in the energy range 31-39 meV. 

Discommensurations could also be responsible for the 
substantial widths of the spin wave modes. In a region of 
commensurate checkerboard order a spin wave can prop- 
agate freely, but when it encounters a discommensura- 
tion the uniformity of the magnetic order is interrupted 
sharply, which could scatter the spin wave. In this case, 
the correlation length of the spin wave peaks should cor- 
respond roughly to the average spacing between discom- 
mensurations. Earlier we found the correlation length to 
be 24 ± 2 A for energies in the range 25-30 meV, which 
corresponds to approximately 9et/\/2, i.e. 9 Ni positions 
when projected along the diagonal of the square lat- 
tice, which is the distance between discommensurations 
according to the model proposed by Kajimoto et al.7 
Further support for the broadening mechanism proposed 
here can be found from a comparison of the widths of 
spin wave peaks observed in neutron scattering measure- 
ments of La2^ x Sr x Ni04 with x — 0.275 and x = 1/3. 9 
The x = 1/3 spin wave peaks show no measurable broad- 



9 



ening, whereas the x — 0.275 peaks are broadened. This 
is consistent with our broadening mechanism since the 
stripe order of x = 1/3 is commensurate and without 
discommensurations, whereas that of x = 0.275 is in- 
commensurate and does have discommensurations. 

This leaves us one remaining feature of the spin exci- 
tation spectrum to consider, namely the low energy dif- 
fuse scattering distributed around the magnetic order- 
ing wavevectors, shown in Fig. As mentioned ear- 
lier, this scattering has a different temperature depen- 
dence to the spin wave scattering, and from this we can 
conclude that it arises from a different magnetic compo- 
nent to the ordered spins. Low energy diffuse scatter- 
ing was also observed in the spin excitation spectrum of 
La5/3Sr 1/ / 3 Ni04»i^ In that case the diffuse scattering was 
almost one-dimensional and, like the present case, there 
was no static component. 

We attributed the quasi- ID diffuse scattering from 
La 5 / 3 Sr 1 / 3 Ni04 to dynamic AFM correlations among the 
S = \ spins attached to lines of Ni 3+ holes in the 
charge stripes, and the origin of the diffuse scattering 
from La 3 / 2 Sr!/ 2 Ni04 is probably the same. Because the 
Ni 3+ holes are arranged on a near-checkerboard pattern 
in La 3/ / 2 Sr 1 / 2 Ni04 the spin correlations among the Ni 3+ 
sites are expected to be quasi-2D, consistent with the 
observed diffuse scattering. The slight elongation of the 
scattering distribution is consistent with the introduc- 
tion of a stripe-like texture into the checkerboard pattern 
by discommensurations. The width of the peak shown 
in Fig. Cfc) is roughly twice that of an equivalent cut 
through the diffuse scattering from La 5 / 3 Sr 1 / 3 Ni04. This 
indicates that the correlations between the Ni 3+ spins are 
weaker in La 3 / 2 Sr]y 2 Ni04. The interpretation of the dif- 
fuse scattering in La 3 / 2 Sriy 2 Ni04 presented here implies 
that there exists an AFM coupling between the Ni 3+ 
spins that generates short-range fluctuations towards a 
checkerboard ordering. The absence of static magnetic 
order on the Ni 3+ sites could have implications for re- 
cent predictions of orbital ordering in La 3 / 2 Sr!/ 2 Ni04^i 

Finally, it is worth comparing our results with 
those obtained on other checkerboard charge-ordered 
compounds. To our knowledge, the only other such 
compound whose spin excitation spectrum has been 
measured in detail is La 3 / 2 Sr 1 / 2 Co04iSS The spin-charge 
order is much closer to a perfect checkerboard pattern 



in the half-doped cobaltate than in the half-doped 
nickelate. The incommensurability observed at low 
temperatures for La 3 / 2 Sr 1 / 2 Co04 is e = 0.49)* compared 
with e = 0.44 for La 3 / 2 Sr!/ 2 Ni04. Recall that for a 
perfect checkerboard order e = 0.5. The measured 
spin excitation spectrum of La 3 / 2 Sr 1 / 2 Co04 exhibits 
a simple spin wave dispersion extending up 16 meV, 
and a second mode at energies around 30 meV which 
is relatively flati^ The spin wave dispersion does 
not exhibit any of the unusual features found in the 
case of La 3 / 2 Sr 1 / 2 Ni04. This can be understood if, 
as argued above, the deviations from a simple spin 
wave picture are due to discommensurations: the very 
small incommensurability of the cobaltate implies very 
few discommensurations are present. The absence of 
low energy diffuse scattering from the cobaltate also 
makes sense since the Co 3+ ions are believed to have a 
non-magnetic singlet ground stated 



V. CONCLUSIONS 

The spin excitation spectrum of La 3 / 2 Sr!/ 2 Ni04 has 
been found to contain a number of interesting features. 
We have argued that the low energy diffuse scattering, 
which resembles a similar signal previously observed in 
measurements on La5/ 3 Sr!/ 3 Ni04, originates from anti- 
ferromagnetic correlations among spins attached to the 
Ni 3+ ions. We have also argued that other strange fea- 
tures in the excitation spectrum of La 3 / 2 Sr 1 / 2 Ni04, such 
as the probable magnetic modes of scattering dispersing 
from (0.5, 0.5) and (1, 0) type positions and the large in- 
trinsic widths of the spin excitations, can be understood 
(at least at the level of hand waving) in terms of a dis- 
commensuration model. The main obstacle in the way of 
a more quantitative account of the spin excitation spec- 
trum of La 3 / 2 Sr 1 / 2 Ni04 is the lack of a detailed model 
for the ground state order. 
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